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Abstract. Debris disks are believed to be related to planetesimals left over around stars after planet formation has ceased. The 
frequency of debris disks around M-dwarfs which account for 70% of the stars in the Galaxy is unknown while constrains have 
already been found for A- to K-type stars. We have searched for cold debris disks around 32 field M-dwarfs by conducting 
observations ai A = 850pm with the SCUBA bolometer array camera at the JCMT and ai A = 1.2mm with the MAMBO array 
at the IRAM 30-m telescopes. This is the first survey of a large sample of M-dwarfs conducted to provide statistical constraints 
on debris disks around this type of stars. We have detected a new debris disk around the M0.5 dwarf GJ842.2 at A = 850yum, 
providing evidence for cold dust at large distance from this star (~ 300 AU). By combining the results of our survey with the 
ones of Liu et al. (2004), we estimate for the first time the detection rate of cold debris disks around field M-dwarfs with ages 
between 20 and 200 Myr. This detection rate is 13^8% and is consistent with the detection rate of cold debris disks (9 - 23%) 
around A- to K-type main sequence stars of the same age. This is an indication that cold disks may be equally prevalent across 
stellar spectral types. 

Key words. Stars : circumstellar matter; surveys; stars: low-mass; planetary systems : formation 
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dwarfs have planets with masses smaller than Jupiter. However 
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revealed substructures that are best interpreted as dust trapped 
in resonance with an unseen planet in the inner part of the sys- 
tem (Wyatt 2003, 2006). Also the large reflecting surface area 
of the dust grains is responsible for scattered light that provides 
a wealth of information useful to disclose the structure and dy- 
namics of these systems (e.g. /3 Pic in Golimowski et al. 2006). 
Cold dust in the Kuiper Belt and warm dust in the asteroid belt 
are the debris of our solar system. 

The first attempt to search for mid-IR excesses from field 
M-dwarfs, i.e. M-stars outside of star forming regions, yielded 
only three detections in the IRAS point source catalogue 
(Inseok Song et al., 2002). However, debris disks around M- 
dwarfs may be missed by such mid-IR surveys because M- 
dwarfs are underluminous (L = 0.1 - 0.001 Lg) so that the 
irradiated dust is cold (<20 K) and more easily detected at 
(sub)mm wavelengths. In this paper, we report on observations 
of 32 young M-dwarfs sA A = 850jum with the Submillimeter 
Common User Bolometer Array (SCUBA) at the James Clerk 
Maxwell Telescope (JCMT), and at i = 1 .2mm with the Max- 
Planck Millimeter Array (MAMBO) at the IRAM 30m tele- 
scope. Previously, only three field M-dwarfs have been ob- 
served in the (sub)mm (Liu et al. 2004). Our aim is to determine 
the frequency of debris disks around M-dwarfs. 

Several studies have adressed the problem of planet for- 
mation around M-dwarfs and their associated debris disks. 
Mechanisms that might be specific to this stellar spectral type 
have been discussed by Johnstone et al. (1998), Laughlin et 
al. (2004), Throop & Bally (2005), Plavchan, Jura & Lipscy 
(2005). Our contribution to establishing the statistics of debris 
disks around M-dwarfs provides new observational constrains. 

Section 2 briefly describes the target stars. Section 3 de- 
scribes the observations at the JCMT and IRAM30m (sub)mm 
telescopes. Section 4 describes the results of these observations 
and provides details on the discovery of the debris disk around 
the MO. 5 dwarf GJ842.2. We draw statistical conclusions on 
debris disks around M-dwarfs in section 5. 

2. Target stars 

Observations have shown that debris disks around main se- 
quence stars of spectral types A to K with ages less than 150 - 
400 Myr are dustier than older ones (Habing et al. 2001, Rieke 
et al. 2005). We therefore selected the youngest M-dwarfs in 
their post planetary formation phase with ages older than 10 
Myr. We shall not denote them main-sequence stars because an 
MO star, e.g., reaches the main sequence only after 900 Myr 
(Siess, Dufour & Forestini 2000). We selected the M-dwarfs 
that belong to the Moving Groups identified by Montes et al. 
(2001) and Zuckerman & Inseok Song (2004a,b), since mem- 
bership to such a group is a reliable age criterium. The se- 
lected stars belong to the Local Association open cluster (age 
= 20 - 150 Myr), IC2391 (35-55 Myr), the AB Dor moving 
group (100 - 125 Myr as revised by Luhman et al. 2005), the 
Castor moving group (200 Myr), the Ursa Maj moving group 
(500 ±100 Myr as revised by King et al. 2003), and the Hyades 
open cluster (600 Myr). Note that HIPl 14066 is part of the AB 
Dor group according to Zuckerman & Inseok Song (2004a,b) 
but is named GJ9809 and assigned to the Local Association 



by Montes et al. (2001). In our surveys, the completness for 
M-dwarfs in the Local Association is 61 % (11 M-dwarfs ob- 
served/18 M-dwarfs in cluster), in Castor it is 75 % (6/8), in 
Ursa Major it is 40 % (4/10), in AB Dor it is 71 % (5/7), in 
Hyades it is 29 % (5/17), in IC2391 it is 50 % (1/2). The M- 
dwarfs of these clusters that are not included in our surveys 
were inaccessible to the telescopes because of their low decli- 
nations. The selected stars have distances between 2.6 pc and 
33.8 pc. Our surveys at A =1.2mm and 850;um include the ob- 
servations of 32 different stars ; 12 were observed at the JCMT 
and 24 at the IRAM 30m telescope with four stars in conmion 
in the two surveys (GJ285, GJ393, GJ9809, GJ4247). 

3. Observations 

3.1. JCMT/SCUBA observations 

The SCUBA bolometer array (Holland et al. 1999) at the JCMT 
in Hawaii at altitude 4092m was used to observe the 12 M- 
dwarfs Usted in Table 1. SCUBA uses two arrays of 37 and 91 
bolometers to simultaneously observe at 850//m and 450jum, 
respectively, the same region of the sky which is ~ 2.3' in 
size. The arrays are arranged as concentric hexagonal rings of 
bolometers at each wavelength. A fully sampled map of the 
2.3' region can be made by "jiggling" the arrays according to 
an optimized pattern (Holland et al. 1999) but requires a long 
integration time to reach the 2 mJy/14"beam sensitivity sought 
in our program. Standard photometry with a single bolome- 
ter of the array (i.e. on-off observations for sky background 
substraction) yields sensitive observations in less time but the 
target source must be smaller than the beam (FWHM=14" at 
850//m). This is not certain for most of our targets because of 
the large a priori uncertainty on disk sizes (factor 3 - 5). Thus, 
in order to optimize the use of our observing time, SCUBA 
data were taken in the wide photometry mode of JCMT de- 
scribed by Sheret, Dent and Wyatt (2004). It ensures that a 
source that is partially resolved or off'set from the observed po- 
sition is detected with maximum sensitivity. In practice, the ar- 
rays are moved through a 12-point pattern that makes the cen- 
tral bolometer of each array cycles through 4 inner points at 1" 
from the target position and 8 outer points at 7". 

In general, an observation consisted of 78 integrations of 
38 seconds, which corresponds to the total on+ofF time of 50 
minutes. During each integration, the telescope beam cycled 
through the 12-point pattern described above and the secondary 
mirror continuously chopped the sky over 60" at 7.8125 Hz 
in azimuth to remove any large-scale sky variations. Each ob- 
servation of a target was preceded and followed by focusing 
and pointing. Atmospheric opacity T225ghz was measured by 
radiometry several times per minute and by "sky dips" once ev- 
ery few hours ; T225ghz> the zenith opacity at 225 GHz, the oper- 
ating frequency of the radiometer, was between 0.06 and 0.08 
during our observations, corresponding to precipitable water 
vapour colunms of ~ 1.2-1.6 nrai. The absolute flux density 
scale was checked by 1 1 measurements of standard caUbrators 
at 850jum (Uranus, Saturn, CRL2688, CRL618) over the 7 ses- 
sions of our program between August 26 2004 and October 07 
2004. From these measurements, we derived a conversion fac- 
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tor FCF - 222 ± 15 (Jy/Volt) for this period. The uncertainty 
in this factor provides the accuracy of our absolute flux den- 
sity scale, which is better than 10 % at 850yum. Observations 
at 450yum are ~15 times less sensitive than at 850jum while the 
flux density is only 6 times higher for dust with a modified 
black-body spectrum and P - \. None of our targets were de- 
tected at this shorter wavelength. 

The SURF package (Jenness et al. 2002) was used to com- 
bine the data from difl^erent integrations, remove anomalous 
spikes, flatfield the array, subtract the sky-background noise 
level and apply atmospheric extinction corrections. The cali- 
brated data were then analysed by a specific software to av- 
erage the data of the central bolometer and to produce the 
map with the data of all the bolometers. This map is incom- 
plete however because the sky is undersampled in the wide 
photometry mode. Nonetheless it allows to reaUably detect a 
source, although interpretation of its measured flux density and 
morphology is somewhat complicated. For the two sources de- 
tected (GJ842.2 and GJ696), we also re-analysed the data with 
the Edinburgh custom data-reduction software written in IDL 
and found similar results. The observations of the 12 M-dwarfs 
were conducted on 2004 August 26, 29, 3 1, on 2004 September 
19 and on 2004 October 3, 6, 7. In the course of this pro- 
gramme, GJ842.2 was observed 4 times, GJ82, GJ212, GJ696, 
and GJ890 twice and the other stars once. 

3.2. IRAM 30m/MAMBO observations 

The 117-channel MPIfR bolometer array MAMBO-2 (Kreysa 
et al. 1998) at the IRAM 30-m telescope on Pico Veleta at 
2900m altitude near Granada in Spain was used to observe the 
24 M-dwarfs Usted in Table 2. MAMBO-2 operates at an ef- 
fective frequency of 250GHz {A= 1 .20mm) with a half-power 
spectrum bandwidth of 80GHz. At /I = 1 .20 mm, the IRAM 30- 
m telescope has an efi'ective beam of 10.7"(FWHM), smaller 
than the angular sizes of the potential debris disks for most 
of our targets. Simple on-off" photometry pointed at the target 
would therefore not be a good strategy to detect a disk and 
multiple pointed observations, similar to wide photometry at 
JCMT, would be inefficient at the IRAM 30-m telescope be- 
cause of the overhead due to the motion of this large telescope. 
Thus, we mapped the targets. 

We used the standard on-the-fly scanning mode, where the 
telescope scans in the azimuthal direction. The signal from the 
sky was modulated by the secondary mirror wobbling over a 
throw of 60" at 2 Hz in the scan direction (azimuth). The 
wobble frequency reflects a compromise between eliminating 
changes in the atmosphere on as short a time-scale as possi- 
ble and the challenges involved in moving a 2-m secondary 
at this frequency and keeping it mechanically stable. For each 
target, we made a map, 400 x 320 arcsec^ in size, scanned at 
a velocity of 4 arcsec s"' and with an elevation spacing of 4 
arcsec. This results in a fuUy-sampled map over this area with 
a ~2mJy/l l"beam rms noise in the central region per observa- 
tion of 33 minutes (on-i-off times). Some targets were observed 
twice or more in difi'erent scan directions and maps were co- 
added to improve this rms noise. In practice, however, weather 



conditions for these multiple observations were diff'erent (opac- 
ity 0.2 < T225GHZ < 0.45, low and medium sky-background 
noise) and thus, quite often, one of the maps dominates the 
final summation (except for GJ628). The opacity of the atmo- 
sphere was measured every other hour by performing a skydip, 
but was also continuously monitored with a radiometer located 
next to the telescope. Short on-off observations of position cali- 
brators were carried out before and after each observation of the 
programm sources to check the telescope pointing which was 
found to be excellent (~ 1"). In order to set the absolute flux 
density scale, primary flux calibrators (including planets) were 
observed at the beginning and end of each run, which resulted 
in an absolute flux calibration uncertainty of ~ 15% (Icr). 

The IRAM data were reduced using the MOPSIC software 
package (Zylka 1998). Data were despiked, filtered for atmo- 
spheric variations, flat-fielded and corrected for atmospheric 
opacity and turned into intensity maps and signal-to-noise ratio 
maps. The observations of the 24 M-dwarfs listed in Table 2 
were conducted in December 2004, in January, February and 
May 2005, and in January and March 2006. Integration times 
were between 0.5 and 2 hours (see Table 2). GJ628 has been 
a special case, after detecting possibly significant peaks in the 
first maps, we extended our observations to 20 hours of inte- 
gration to reach a map noise rms of 0.4mJy/l l"beam in order 
to search for diffuse emission between these peaks. 

The relevant sensitivities for our SCUBA and MAMBO 
surveys are in terms of circumstellar dust mass rather than mJy. 
Although the sensitivity is ~ 2 mJy/beam for each observations 
in both surveys, the sensitivity with MAMBO was enhanced 
by taking a higher number of observations and by averaging 
the data over a canonical disk area assumed as large as the 
Kuiper-belt. This averaging can be done also with the SCUBA 
wide photometry data but not as efficiently because the sky is 
not fully sampled. Also the dust mass sensitivity depends on 
the steUar distance squared, and the M-dwarfs observed with 
MAMBO are closer than the ones observed by SCUBA in gen- 
eral. The all-in sensitivity is ~ 1 lunar mass of circumstellar 
dust for both surveys. 

4. Results 

4.1. JCMT/SCUBA results 

In our JCMT/SCUBA sample of 12 M-dwarfs, we have 
detected significant peaked and extended emissions around 
GJ842.2 revealing a debris disk around this M0.5 dwarf. We 
analysed the data of GJ842.2 in two ways. First, we straight- 
forwardly averaged the data of the central bolometer and made 
a significant detection with a flux density of 4.8 ±1.3 mJy. In 
efi'ect, this flux density is the average of the emission at the 4 
points 1" away from the star position and at the 8 points 7" 
away and may be underestimated depending on the morphol- 
ogy of the source (Sheret, Dent and Wyatt, 2004). Nonetheless, 
the source is detected and this is the important point for the sur- 
vey. Second, we averaged the data from the first ring of bolome- 
ters (6 bolometers at mean radius 25" from center), and found 
the significant mean brightness of 2.4 ± 0.6 mJy/14"beam 
over this area. We have also averaged the data of the second 
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Table 1. JCMT/SCUBA wide photometry at 850yum. 



Star 


Sn 


RA" 


DEC 


Dist. 


Lo2(L,)* 


Moving 


Integration 


Flux density 


Size " 


Dust 


Dust 




tvne 


(J2000) 


(J2000) 




Lofff erffs/s^ 


Group 


time 






temp'' 


mass^ 
























(moon) 


GJ82 


M4 


01 59 23.5 


58 31 16 


12.0 


28.75 


Loc Ass. 


2 


2.0 ± 1.4 


84 


16 


< 1.3 


GJ212 


M0.5 


05 41 30.7 


53 29 20 


12.5 


28.57 


Loc Ass. 


2 


1.3 ± 1.4 


88 


23 


< 0.8 


GJ285 


M4.5 


07 44 40.1 


03 33 08 


5.9 


28.62 


Loc Ass. 


1 


-0.7 +1.9 


41 


20 


< 0.3 


GJ393 


M2 


10 28 55.5 


00 50 27 


7.2 


26.84 


Loc Ass. 


1 


1.9 ± 1.9 


50 


26 


< 0.3 


GJ507.1 


Ml. 5 


13 19 40.1 


33 20 47 


17.4 


n/a 


Loc Ass. 


1 


-0.4 ± 2.0 


121 


19 


< 2.8 


GJ696 


MO 


17 50 34.0 


-06 03 01 


21.9 


n/a 


Loc Ass. 


2 


0.8 ±0.8 


150 


20 


< 1.9 


GJ9809 


MO 


23 06 04.8 


63 55 34 


24.9 


29.43 


Loc Ass. 


1 


-5.2 ± 2.3 


174 


18 


<7.6 


GJ4247 


M4 


22 01 13.2 


28 18 25 


9.0 


n/a 


Castor 


1 


1.1 ±2.1 


63 


18 


<0.9 


GJ277B 


M3.5 


07 31 57.3 


36 13 47 


11.5 


29.05 


Castor 


1 


-2.1 ± 1.8 


80 


17 


< 1.3 


GJ842.2 


M0.5 


21 58 24.8 


75 35 21 


20.9 


n/a 


Castor 


4 


25 ± 4.6^ 


300 


13 


28 ±5 


GJ890 


M2 


23 08 19.5 


-15 24 35 


21.9 


29.22 


Castor 


2 


-2.6 ± 1.6 


153 


16 


<4.5 


GJ875.1 


M3 


22 51 53.7 


31 45 15 


14.2 


28.92 


IC2391 


1 


0.52 ± 2.1 


99 


17 


<2.4 



Observed positions, i.e. coordinates updated with the Hipparcos proper motions at the mean observation date of September 2004. 
* ROSAT All Sky Survey of Nearby Stars : Huensch et al. 1999. 
' Size is the radius of a debris disk as large as the telescope beam (14"). 

Temperature is computed at this radius. 
' Mass upper Unrits are based on the 3cr flux densities (see text). 
f Total flux density of the extended emission and the South-West peak. 



ring of bolometers further away at 45" and found 0.83 ± 0.70 
mJy/beam, i.e., consistent with no detection. The third ring at 
65" was used to remove sky-background noise and, thus, de- 
fines the zero-level of the flux density scale. This implies that 
extended emission exists out to ~ 25" and so the emitting disk 
is as large as ~ 500 AU at the distance of GJ842.2 (20.9 pc). In 
an attempt to identify the general morphology of this extended 
emission, we split the field of view in eight sectors to average 
the data in each one of them (Fig. 1). In this decomposition, we 
used only two free parameters, the global orientation of the sec- 
tors and the arclength of the inner South-West sector, maximiz- 
ing the signal-to-noise ratio in it. The sector radii are not free 
parameters but imposed by the SCUBA array. This averaging of 
the data indicates that the emission is non-uniformly distributed 
and is oriented along the North East - South West direction, 
suggesting an incUned disk. The plot at the bottom of Fig. 1 is a 
section of the brightness along this direction. The total flux den- 
sity (brightness X surface area) is estimated to be ~ 25±4.6mJy 
(i.e. S/N=5.5) by adding the flux densities of the 4 sectors along 
the Norh East - South West direction. Finally, we computed the 
hkehhood that the offset peak corresponding to the South-West 
sector (~ 11 mJy at 7") is a background source. The probabil- 
ity to find a background source within a radius of 7" from an 
observed position in our 12 star survey is as low as 3%, com- 
puted from 2100 X (1 l/2.3)-i-5 X na"/3600"f x 12, based on 
the source count power-law N(S) oc S^^'^ normalized by 2100 
sources/deg^ for S > 2.3 mJy at 850//m (Rowan-Robinson 
2001). So we conclude that there is a cold debris disk around 
GJ842.2, possibly inclined, made of a peak of emission at 7" 
Soufli West from the star (~ 140 AU at 20.9pc) embedded in 
low brightness emission extending over 25" (~ 500 AU). We 
emphasize that this attempt to identify the general morphology 



of the disk is speculative and that the straightforward averag- 
ing of the data of the first ring of bolometers mentioned above 
that indicates an extended source of mean brightness 2.4 ± 0.6 
mJy/14"beam is our main result. 

The total SCUBA flux density Ssso^m of GJ842.2 was 
used to calculate the dust mass Mj of the disk by assum- 
ing a simple model of large grains (lOOjum in size) located 
at the mean radius of the extended disk (300 AU). We use 
the standard optically thin dust formula for (sub)mm emission 
\Sa - Md X B{A,Tg) X Kabs I d^, where Tg is the grain 

temperature, d is the star distance, Kabs = Ksso^m ( ^^"^"' ) 
is the mass opacity for the modified black body, with 
KssOfim = 1-7 cm^ g"' and j8 = 0.8 for ~ lOO/im size 
grains. Dent et al. (2000)]. At 300 AU, the dust temperafiire 
is only 13 K and the total dust mass is 28 + 5 lunar masses for 
Sgsoimi = 25 ±4.6 mJy. Those simple assumptions provide only 
the magnitude of the cold dust mass for the disk of GJ842.2 but 
it shows clearly that it is larger than the dust mass estimated 
around the Ml dwarf AU Mic (~ 1 lunar mass) and around the 
M0.5 dwarf GL182 (~ 2.1 lunar mass) by Liu et al. (2004). 

In addition, we reanalysed IRAS data of GJ842.2 with the 
on-line program scanpi and determined the color-corrected flux 
density 57± 12 mJy at 12//m. This flux density matches the pho- 
tospheric emission computed from the NextGen model for this 
MO dwarf (Tphoto = 3500 K, g = 5.5, [Fe/H] = 0.0) by Allard et 
al. (2000, 2001), normahzed by the B,V,R,I,J,H,K flux densities 
of GJ842.2 accessible in the SIMBAD data base (B, V bands : 
Third Catalogue of Nearby Stars ( 1 99 1 ) by Gliese and Jahreiss, 
Astron. Rechen-Institut, Heidelberg ; B, R, I bands : USNO- 
Bl.O catalogue, Monet et al. 2003 ; J, H, K bands : 2MASS 
All-Sky catalogue of point sources, Skrutskie et al. 2006). The 
scanpi signals in the other IRAS bands were not significant ex- 



Lestrade et al.: Search for Debris Disks around M-dwarfs 



5 




50 , ^ -50 

distance from star (arcsec) 



Fig. 1. Brightness distribution around the MO. 5 dwarf GJ842.2 
atA = 850/im based on the data from all the bolometers of the 
JCMT/SCUBA array. We made an attempt to identify the gen- 
eral morphology of the source in splitting the field of view in 
eight sectors and averaging the data in each one of them. As 
it is seen, the emission is non-uniformly distributed and is ori- 
ented along the North East - South West direction, suggesting 
an incUned disk. The plot at the bottom of the figure is a sec- 
tion of the brightness along this direction. The total flux density 
(brightness x surface area) is estimated to be ~ 25+4.6mJy (i.e. 
S/N=5.5) by adding the flux densities of the 4 quadrants along 
this direction. The FWHM beamwidth of the telescope is 14" 
at Observed position of GJ842.2 is a = 21h58m24.8s 

and 5 = 75°35'21" (J2000). 

cept possibly at 25/im (58 ±13 mJy, color corrected). However, 
an IRAS flux density must be treated circumspectly when the 
signal-to-noise ratio is lower than 5. If this 25 flux density 
(4.5 <t) were real however, it would be clearly higher than the 
photospheric level of GJ842.2 and so, interestingly, indicative 
of warm dust close to the star. This needs confirmation by new 
observations. 

For the MO. 5 dwarf GJ696 in our survey, the central 
bolometer data did not yield any significant flux density in the 
central region close to this star (radius < 14") and data aver- 
aging of the first ring of bolometers (18" - 35") did not reveal 
any extended emission at the Icr level of 0.8 mJy/14"beam. 
However, we have found a significant peak of emission, 1 1 mJy 
(4o-) at 20.8" and PA=-h73° (East-Norfli-East) from the posi- 
tion of this star in the SCUBA incomplete map. It it is possible 
that its true flux density is underestimated because of under- 



sampling as mentioned above. This peak is 10 times higher than 
the noise of the first ring of bolometers. Thus, we do not regard 
this detection as evidence for a debris disk around GJ696, re- 
ferring to the SCUBA map of e Eri which shows ratio between 
peaks and difl'use emission of less than 2 (Greaves et al. 2005). 
It is likely that the peak near GJ696 is a background source. 
We computed the likelihood a sub-mm background source be 
within a radius of 21" over the whole 12 star sample using the 
source count at 850jum by Rowan-Robinson (2001) and found 
a probability of 25 %. There is no corresponding 1 .4GHz radio 
counterpart in the NVSS catalogue by Condon et al. (1998) and 
no optical object in the USNO-B 1 catalogue that is positionally 
coincident. 

No peaked or diffuse emissions have been detected around 
the other 10 stars of our survey. In Table 1, we derived the upper 
limits on the dust masses of their potential debris disks using 
the 3cr flux density upper Umits and the dust temperature for a 
linear separation from the star corresponding to 7" which is the 
angular extent well probed by wide photometry. 



4.2. IRAM 30m/MAMBO results 

In our IRAM/MAMBO sample of 24 M-dwarfs (Table 2), we 
have detected emission in form of 5 isolated sources around 
the M3.5 dwarf GJ628 (Fig. 2). We have investigated whether 
or not these sources might simply be background galaxies. We 
have used the large catalogue USNO Bl.O (Monet et al. 2003) 
containing more that a billion objects with position uncertain- 
ties of ~ 0.3" to search for optical counterparts and sum- 
marized our results in Table 3. The only convincing associ- 
ation is for the source named MAMBO-W at ~ 154" West 
of GJ628 in Fig. 2, based on coincidence between MAMBO 
and USNO positions at the Icr level. In addition, this source 
MAMBO-W matches within ~ 1 .Icr the position of the radio 
source NVSS163007-123940 in the 1.4 GHz NVSS catalogue 
by Condon et al. (1998), if cr is the NVSS and MAMBO posi- 
tion uncertainties quadratically combined. Hence, we conclude 
that MAMBO-W is a background source identified in the vis- 
ible and in the radio as a quasi-flat spectrum source from cm 
(NVSS : 5.6 ± 0.5 mJy) to millimeter wavelengths (this pa- 
per : ~ 3.6 mJy). For the other four MAMBO sources, the 
likelihood of optical associations is low (position discrepan- 
cies in Table 3) and they have no radio counterpart in the NVSS 
catalogue. Using the density of mm background sources, 600 
sources deg"^ for Sa > 2 mJy at 1 .2 mm, found by Voss et al. 
(2006), we expect 1.1 background sources over the field that 
includes these four sources (r < 90", dotted circle in Fig 2). 
In this probability calculation, we did not include the fields of 
the other 23 stars observed because their maps are not as deep 
as the one of GJ628. Since we observe more MAMBO sources 
than statistically expected from the sky background, they may 
be instead four dust clumps that are part of a large disk (ra- 
dius ~ 400 AU) surrounding the M3.5 dwarf GJ628. In this 
case, their temperatures would be between 8 and 12K and their 
dust masses would be ~ 0.1 lunar mass each. We found no sig- 
nificant mm diffuse emission between these point-like sources 
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at the 0.4 mJy/1 l"beam level after 20 hours of integration. If 
GJ628 were similar to e Eri, the expected diffuse component 



Table 2. IRAM/MAMBO aperture photometry at ^ = 1.2mm. 



Star 


Sp. 


RA" 


DEC 


Dist. 


LogCL,)* 


Moving 


Integr. 


e 


Mean brightness " 


Map noise rms 




Dust 


Dust 




Type 


(J2000) 


(J2000) 


(pc) 


Log(ergs/s) 


Group 


time 
(hrs) 


(") 


(mJy/ll"beam) 


(mJy/ll"beam) 


(30-) 
(mJy) 


temp^ 
(K) 


mass f 
(moon) 


GJ65A 


M5.5 


01 39 02.6 


-17 56 59 


2.6 


27.59 


Hyades 


1.1 


46 


-0.30 ±0.5 


2.1 


26 


15 


<0.9 


GJ109 


M3 


02 44 15.8 


25 31 22 


7.6 


27.30 


Hyades 


1.3 


16 


0.67 ±1.4 


2.1 


9 


21 


< 1.6 


HIP16563B« 


MO 


03 33 14.0 


46 15 19 


33.8 


n/a 


AB Dor 


0.9 


4 


0.35 ±1.0'' 


- 


3 


27 


< 8.0 


HIP 17695 


M3 


03 47 23.4 


-01 58 21 


16.3 


n/a 


AB Dor 


0.5 


5 


-0.79 ±2.2 


2.2 


7 


20 


<5.6 


GJ3379 


M4 


06 00 04.5 


02 42 23 


5.4 


n/a 


Hyades 


0.5 


22 


-0.06 ±1.1 


2.3 


14 


18 


< 1.5 


GJ234A 


M4 


06 29 23.6 


-02 48 53 


4.1 


n/a 


Loc Ass 


1.1 


29 


0.76 ±0.7 


1.8 


14 


18 


<0.9 


GJ285 


M4.5 


07 44 40.0 


03 33 06 


5.9 


28.62 


Loc Ass 


1.1 


20 


-0.13 ±1.0 


1.7 


9 


17 


< 1.1 


GJllll 


M6.5 


08 29 49.1 


26 46 29 


3.6 


n/a 


Castor 


0.5 


33 


-0.34 ±0.8 


2.4 


22 


12 


< 1.9 


GJ393 


M2 


10 28 55.3 


00 50 24 


7.2 


26.84 


Loc Ass 


1.7 


18 


-0.38 ±1.1 


1.5 


7 


24 


<0.8 


GJ402 


M4 


10 50 51.8 


06 48 25 


6.8' 


n/a 


Loc Ass 


1.1 


18 


-0.24 ±1.0 


1.7 


8 


18 


< 1.0 


GJ447 


M4 


1147 44.5 


00 48 10 


3.3 


26.84 


Ursa Ma 


1.1 


36 


0.28 ±0.4 


1.4 


14 


18 


<0.6 


GJ408 


M2.5 


11 00 04.1 


22 49 57 


6.6 


26.60 


Castor 


1.1 


18 


0.16 ±1.0 


1.7 


8 


22 


< 1.0 


GJ569A 


M2 


14 54 29.3 


16 06 03 


9.8 


28.54 


Ursa Ma 


1.1 


12 


1.76 ±2.2 


2.2 


7 


24 


< 1.8 


GJ628 


M3.5 


16 30 18.0 


-12 39 51 


4.3 


26.48 


Loc Ass 


20.0 


28 


0.1 ±0.2 


0.4 


3 


12 


<0.4 


GJ625 


M1.5 


16 25 24.8 


54 18 13 


6.6 


26.70 


Ursa Ma 


0.5 


18 


-0.85 ±1.5 


2.5 


12 


26 


< 1.3 


HIP81084 


M0.5 


16 33 41.5 


-09 33 12 


31.9 


n/a 


AB Dor 


0.5 


4 


0.63 ±2.6 


2.6 


8 


27 


< 19 


HIP86346 


MO 


17 38 39.6 


61 14 16 


24.5 


n/a 


AB Dor 


0.5 


5 


-0.13 ±2.1 


2.1 


6 


27 


< 8.2 


GJ791.2 


M4.5 


20 29 48.6 


09 41 19 


8.9 


27.89 


Hyades 


1.1 


14 


-0.98 ±1.7 


1.7 


6 


17 


<2.0 


GJ849 


M3.5 


22 09 40.7 


-04 38 26 


8.8 


n/a 


Hyades 


0.5 


14 


-1.32 ±2.7 


2.7 


10 


19 


<2.7 


GJ876 


M4 


22 53 17.0 


-14 15 52 


4.7 


n/a 


Loc Ass 


0.5 


25 


-0.07 ±1.5 


2.8 


19 


18 


< 1.6 


GJ873 


M3.5 


22 46 49.3 


44 19 59 


5.0 


29.09 


Ursa Ma 


0.6 


24 


0.32 ±1.1 


2.0 


13 


19 


< 1.1 


GJ4247 


M4 


22 01 13.2 


28 18 25 


9.0 


n/a 


Castor 


1.1 


13 


-0.18 ±1.4 


1.4 


5 


18 


< 1.5 


GJ856A 


M3 


22 23 29.1 


32 27 32 


16.1 


29.46 


AB Dor 


1.1 


8 


-0.56 ±1.5 


1.5 


4 


21 


<3.3 


GJ9809 


MO 


23 06 04.8 


63 55 34 


24.9 


29.43 


Loc Ass 


0.5 


7 


0.36 ±2.2 


2.2 


7 


26 


< 10 



" Observed positions, i.e. coordinates updated with the Hipparcos proper motions at the mean observation date of July 2005. 

* ROSAT All Sky Survey of Nearby Stars : Huensch et al. 1999. 

' Mean brightness computed over the effective area of angular size corresponding to the canonical disk 120 AU in diameter. 

Total flixx density Fg is the 3a- map noise integrated over this canonical disk. 

Dust temperature computed at 60 AU from the central star. 
^ Dust mass upper limits based on F^. 

^ Designation and spectral type from Zuckerman & Inseok Song, 2004b and Lopez-Santiago et al., 2006. 

* On-off single bolometer MAMBO standard photometry. 

' Based on parallax from the Research Consortium on Nearby Stars RECONS project (www.chara.gsu.edu/RECONS). 
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Fig. 2. IRAM/MAMBO Signal-to-Noise map of the field 
around the M3.5 dwarf GJ628 at A = 1.2mm. In this map, 
Signal-to-Noise ratio (S/N) is computed from the intensity 
map reduced with the shift-and-add method using the software 
package mopsic (Zylka 1998) and smoothed to the angular res- 
olution 12". The S/N is computed by dividing the intensity at 
each pixel by the local rms (10 x 10 pixels) which is not uni- 
form across this intensity map made by the scanning technique. 
The local rms cr is ~ 0.4 mJy/1 l"beam in the central region (r 
< 60") and degrades towards the edges of the map ; for ex- 
ample, Icr is 0.75 mJy/ir'beam at r ~ 150". The countours 
are -3cr, -2cr (dotted lines), and 2cr, 3cr, 4cr; The map pixel 
size is 3.5". The West component (MAMBO-W) is associated 
with a 1.4GHz radio source (x symbol) of the NVSS cata- 
logue and with an optical object of the USNO-Bl.l catalogue 
(+ symbol), and so it is a background source. The other four 
MAMBO components (NW, NE, E, SW) may be dust clumps 
of a large debris disk around GJ628 but the lack of diffuse emis- 
sion between them requires additional observations for confir- 
mation. Observed position of GJ628 at center of the map is a = 
16h30ml8.0sand5 = -12°39'51" (J2000). 

should be seen at the 2cr level in Fig. 2. Note that close inspec- 
tion of this map does show that there is a l-cr to 2-cr emission 
extending between the E and NE MAMBO sources. It was not 
possible to image deeper with the MAMBO array in a reason- 
able oberving time to confirm such an extended emission. Only 
additional observations, in the Far-IR, in scattered light or by a 
long-term astrometric monitoring of the mm sources that would 
share the proper motion of GJ628 (~ 1 .2"/yr) if they are part of 
its debris disk, could provide decisive evidence for this intrigu- 
ing result. We don't retain GJ628 as having a disk at this stage 
of our investigation. 

We found the significant IRAS/scanpi flux densities 385 ± 
21 mJy at 12/zm and 110 ± 25 mJy at 25yum (color-corrected) 
for GJ628, but they match the NextGen photospheric emission 
computed for this M3.5 dwarf. So no mid-IR emission excess 
above photospheric level is observed. 

For the other stars, no emission is apparent in their 
MAMBO maps, and no significant mean brightness is detected 



when averaging the data over an efi'ective disk with an angu- 
lar diameter (6) corresponding to the canonical diameter 120 
AU, i.e. the size of the Kuiper Belt. Hence, we have used the 
map noise rms of each star to derive the 3-cr upper limit of its 
dust mass by averaging data over this canonical disk. In such 
a derivation, the 3cr total flux density over the effective area 
(e(")) is = 3 X rms x (6'(")/ll") (mJy), which takes into 
account that the measured brightness uncertainty decreases as 
rms/(0/ir'), i.e. as rms/ Vnumber of beams in the disk area, 
while the total flux density increases as (O/ll")^, i.e. as 
(number of beams in the disk area)^. The dust mass is com- 
puted from this total flux density Fq using the optically thin 
dust model described above substituting Fg for and taking 
the absorption coefficient Kabs = 1.3 cm^g"^ at A = 1.2 mm, 
according to the modified black-body law and /3 = 0.8. The im- 
plicit assumption of a disk seen face-on in this derivation pro- 
vides a conservative upper limit for the total flux density and 
hence for the circumstellar dust mass since an inclined disk 
would be covered by a smaller number of beams. In Table 2, 
we report the values for 0("), mean brightness, brightness un- 
certainty, map noise rms, Fo and dust mass upper limit. We do 
not confirm the 2cr fiux density of 17 mJy for GJ873 (EV Lac) 
atA= 1.1 nam reported by Doyle and Mathioudakis (1991). 

5. Discussion 

5.1. Debris disks in our sample of 32 M-dwarfs 

Our surveys of 32 young M-dwarfs conducted at A = 850;um 
and 1 .2mm yielded the detection of one disk around the MO. 5 
dwarf GJ842.2. The emission from this disk is a peak at 7" 
Soufli West from the star (~ 140 AU at 20.9pc) embedded in 
low brightness emission extending over 25" (~ 500 AU). We 
fitted a simple dust model by assuming large grains (~ lOOjum) 
at mean disk radius (300 AU, 13K) that yielded 28 + 5 lunar 
masses of cold dust for the total flux density measured. No 
other map shows convincing evidence for a debris disk in our 
survey, although the SCUBA incomplete map of the MO dwarf 
GJ696 (21.9 pc) shows a locaUzed source (11 ± 2.7 mJy) at 
21" from the star but no extended emission at the level of 0.8 
mJy/14"beam over 35", and the MAMBO map of the M3.5 
dwarf GJ628 (4.3 pc) shows four mm sources within 90" (400 
AU) from this star but no significant diffuse emission at the 0.4 
mJy/1 r'beam level. We do not consider these isolated clumps 
as sufficient evidence for debris disks around these last two 
stars. However, one of the most defining features of debris 
disks that have been imaged in the (sub)mm is that they are 
not smooth, but clumpy, and so we plan other observations to 
investigate these cases further. 

It is interesting to discuss the lifetimes of grains in the de- 
bris disk of GJ842.2. Due to the low emission efficiency in the 
sub-mm of micron-sized grains, the A = 850/im emission is ex- 
pected to arise from grains larger than lOOyum (Wyatt & Dent 
2002). Their coUisional lifetime is tcoii > 250 Myr and their 
Poynting-Robertson (P-R) drag lifetime > 22 Gyr at distance 
r > 140 AU. These lifetimes are larger than the star age (200 
Myr) and so it means that these large grains at r > 140 AU 
could be original ("primordial") from the early protoplanetary 
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Table 3. Positions of the five IRAM/MAMBO sources in the map of GJ628 at A = 1.2mm (Fig. 2) and their nearest optical 
objects referenced in the USNO-Bl.O catalogue. 





MAMBO 


USNO-Bl.O 


Name 




A" 






WUJCLL-MJUICC 






(J2000) 


(J2000) 


density 


separation 


separation 


magnitudes 








(mJy) 


(") (AU) 


(") 




MAMBO-W 


16 30 07.5 


-12 39 42 


3.6 


154 660 


3.0 (l.Oo-) 


B=20.03 R=19.12 


MAMBO-NW 


16 30 15.6 


-12 39 11 


2.7 


53 230 


13.3 (4.0O-) 


B=20.50 R= 17.90 


MAMBO-SW 


16 30 13.9 


-12 39 59 


2.0 


90 390 


10.1 (3.0O-) 


B=19.27R=18.23 


MAMBO-NE 


16 30 19.7 


-12 39 32 


2.4 


83 360 


9.1 (2.5a-) 


B=20.06 R=18.30 


MAMBO-E 


16 30 22.3 


-12 39 25 


1.8 


68 290 


6.9 (2.0a-) 


B=21.63 1=18.61 



" The coordinates uncertainties are 3.5" for MAMBO. 
*■ Flux densities are the peak values unaffected by gaussian smoothing of the map. 



phase. On the other hand, grains closer to the star than 140 AU, 
if any, could not be primordial because of their shorter tcoii, i-e. 
more frequent constructive or destructive collisions. The large 
disk around GJ842.2 is thus unique in terms of the possible 
presence of primordial and processed large dust grains at such 
a late stage of evolution and impHes a large protoplanetary disk 
at birth. Previous to our survey, only the field M-dwarfs AU 
Mic and GJ182 were detected in the sub-mm (Liu et al. 2004). 
These two stars are 10-20 Myr old (see Favata et al. 1998, and 
Barrado y Navascues et al. 1999 for the age of GJ182). The 
disk of AU Mic was resolved by several workers (Kalas et al. 
(2004), Liu (2004) and Krist et al. (2005)) that found the size ~ 
200 AU smaller than the one for GJ842.2 while these two stars 
have the same spectral type but one is twenty times younger. 
Further investigations should attempt to estabUsh whether the 
environmental factor or the evolutionary factor controls these 
two disks. 

Plavchan, Jura & Lipsky (2005) have not detected any 
warm debris disks via 11/im excess in their sample of 9 M- 
dwarfs older than lOMyr. These authors highlight that stellar 
wind drag rather than P-R drag is responsible for dispersal of 
the dust in these systems and this might explain the dearth of 
debris disks around M-dwarfs. For GJ842.2 (L» = 0.063Lo), 
a stellar wind with a mass-loss rate Mj„ 20 times the solar 
rate (M© ~ 2 10"'"* Mo/yr) would be sufficient to remove the 
primordial grains of its disk in 200 Myr according to the P-R 
to stellar wind drag timescales ratio tpRjtsw °^ Msw/L* (equa- 
tion (3) of Plavchan, Jura & Lipsky (2005)). In fact, winds with 
mass-loss rates 1000 times the solar rate are observed for M- 
dwarfs, especially when young, and would make the existence 
of primordial dust after 200 Myr inconceivable. Unfortunately, 
no X-ray or Call data are available that would allow an esti- 
mate of Msw for GJ842.2. However, GJ842.2 is not a flare star, 
and its Hq, absorption (equivalent width EWh„ = - 0.545 A, 
in Gizis, Reid & Hawley, 2002) indicates a moderate chromo- 
sphere and magnetic activity, and thus only a modest steUar 
wind (MuUan et al. 1992, 2001, and Wargelin & Drake 2001). 



Hence, one can speculate that GJ842.2 does not have a mass- 
loss rate larger than 20 times the solar rate and its distant cold 
grains are unperturbed at age 200Myr. Although GJ628 is not 
retained as having a debris disk in the present analysis, it is in- 
teresting to note that its X-ray luminosity (10^^'^^ ergs s"' mea- 
sured by ROSAT, Schmitt & Liefke (2004)) is low, similar to 
that of the quiet Sun (10^^-^ ergs s"'), and its EWh^ (-0.234 A, 
in Gizis, Reid & Hawley 2002) indicates a weak to moderate 
chromosphere. Hence, similarly to GJ842.2, if GJ628 were sur- 
rounded by cold grains, one can speculate that they would be 
unperturbed at age 20-150Myr. Finally, we note also that there 
is no difl^erence in disk detection rates between the subsample 
of 11 stars with high Lx (i.e. Log(Lj:) > 27.5) and the subsam- 
ple of 7 stars with low in Tables 1 and 2. 



We have critically reviewed the literature to list the proper- 
ties of the debris disks around A- to K-type main sequence stars 
in Table 4 and plotted their dust masses as a function of spec- 
tral types in Fig 3 after adding our results for the M-dwarfs. The 
200 Myr old debris disk around GJ842.2 is relatively massive 
(28 lunar mass of dust) and may be the tip of the iceberg for 
the M-dwarf disks. We note that this mass is comparable to the 
highest disk mass of the A- to K-type main sequence stars plot- 
ted in Fig. 3. Future surveys will have to establish whether or 
not this upper envelope for the non-detected M-dwarfs in Fig. 3 
is a real envelope to assert whether or not disks around M- 
dwarfs are as massive as around A- to K-type stars. Presently, 
this comparison is also complicated by the fact that the debris 
disks detected in the (sub)mm around these A- to K-type stars 
are from an IRAS-selected sample from essentially an all-sky 
survey, while the M-dwarfs in our analysis are from an unbi- 
ased (except for age) sample of 32 stars. EventuaUy, a broader 
comparison with protoplanetary disk masses and central star 
masses (Natta, Grinin & Mannings (2000) and Andrews & 
WiUiams, 2005) will allow to discuss how the mass of the cen- 
tral star impacts planet formation. 
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Table 4. Dust masses of debris disks around A- to K-type main sequence stars. Note that HD34700 and HD39944 previously 
classified as debris disks have been rejected; HD34700 is a T Tauri (Sterzik et al. 2005) and has gas emission (Dent et al. 2005) 
and HD39944 was confused with a background galaxy (Moor et al. 2006). 
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0.85 


60 


64 
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" Greaves et al. 2004 ; * Greaves et al. 1998 ; ' Liu et al. 2004 ; '' Najita & Williams 2005 ; ' Zuckerman et al. 1993 ; Sheret et al. 2004 ; 
* Wyatt et al. (2003) ; ' Williams et al. 2004 ; ' Holmes et al. 2003 ; ' Sylvester et al. 1996 ; ' Wilner et al. 2003 ; Wyatt et al. 2005 



5.2. Detection rate of dusty debris disks around 
M-dwarfs 

We have combined our surveys (32 field M-dwarfs observed 
and 1 detection) with the JCMT survey of 3 field M-dwarfs 
with 2 detections by Liu et al. (2004). Note that their small 
sample has the same characteristics in ages, distances and sen- 
sitivity as ours ; AU Mic (9.9 pc) and GJ799 (10.3 pc) of the 
Pic moving group (age 10 - 20 Myr) and GJ182 (26 pc) 
of the Local Association. In this combination, there are (1+2) 
detections out of (32-h3) M-dwarfs searched, yielding the de- 
tection rate / = for cold debris disks in this sample 
of M-dwarfs with ages between 20 and 600 Myr observed in 
the (sub)mm. We have derived the uncertainty by computing 
the full width half maximum of the probability distribution 
P(k/N) = C*(l -/)'^"*/* that k detections be found in a sample 
of A' stars searched if the a priori frequency of occurence is /. 

In order to compare detection rates between surveys, it 
is important to compare observations in the same wavelength 
regime since mid-IR surveys probe warm dust (~ lOOK) close 
to the star (~ I AU) while (sub)mm surveys probe cold dust 
(~20K) distant from the star (lOO's AU). In the Uterature, there 
are two surveys of A- to K-type stars conducted in the sub-mm 
that are suited for comparison with our detection rate. A sam- 
ple of A- and F-type main sequence stars not detected by IRAS 
with ages 10 - 170 Myr was conducted at 850 pim by Wyatt, 
Dent & Greaves (2003) and a sample of F- to K- type stars 



with ages 10 - 180 Myr was conducted at 850 fjm by Najita 

and Williams (2005). It is important also to trim the samples 
so that they have same age interval since the decay of dust in 
debris disks with time seen for A-type stars (Rieke et al. 2005) 
is likely to affect M-dwarfs too. For these considerations, we 
therefore have restricted our (sub)mm sample to the (3-h20) M- 
dwarfs with ages 20-200 Myr, including the three M-dwarfs 
studied by Liu et al. (2004) satisfying this age limit. We then 
find that our detection rate / becomes 13;|;g% (3/23) for cold 
debris disks surrounding field M-dwarfs with ages between 20 
and 200 Myr. This detection rate / compares well with the de- 
tection rate (9 - 14%) found in the A- and F-type star sam- 
ple of Wyatt, Dent & Greaves (2003), and to the detection rate 
23^1^ % (3 detections for 13 stars observed) in the F- to K- 
type star sample of Najita and Williams (2005). It is remark- 
able that these cold dust detection rates for M-dwarfs and A- to 
K-type dwarfs match. This may indicate that the populations of 
cold debris disks among early and late-type stars are alike. We 
have also considered the sub-sample of 89 F- to K-type stars 
with ages between 10 and 300 Myr observed at A = 1.2mm, 
2.7mm and 3mm by Carpenter et al. (2005) that has a detec- 
tion rate of 2%. It is somewhat lower than the rates above but 
the longer wavelengths used in their observations would have 
required deeper observations to be directly comparable to the 
other surveys. 
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Fig. 3. Dust masses of debris disks (filled circles) around A- 
to K-type main sequence stars (Table 4), around the two M- 
dwarfs AU Mic and GL182 from Liu et al. (2004) and around 
the M-dwarf GJ842.2 reported in this paper In addition, we 
plot the dust mass 3cr upper limits of potential disks for the 
non-detections of this paper and for the M-dwarf GJ799 re- 
ported in Liu et al. (2004) (red thin triangles^ SCUBA and 
green thick triangles= MAMBO). The dust mass for the Kuiper 
Belt is shown here to indicate that present observations pro- 
vide dust masses that are only the envelope of a possibly much 
larger population of debris disks. The dashed lines are plotted 
as a guide. 

6. Conclusions 

We have searched for cold debris disks around 32 field M- 
dwarfs at sub-mm and mm wavelengths. We report evidence 
for a new debris disk around the MO. 5 dwarf GJ842.2 of age 
~ 200 Myr which is massive (28 + 5 lunar masses) and large 
(~ 500 AU). We report also the detection of four mJy sources at 
1.2mm around the M3.5 dwarf GJ628 that might be associated 
with a debris disk around this star unless they are background 
sources. 

By combining the result of our surveys with the two cold 
debris disks discovered by Liu et al. (2004) in their observa- 
tions of three field M-dwarfs, we derive for the first time the 
detection rate / = 13^^% for cold debris disks surrounding 
M-dwarfs with ages 20-200 Myr This detection rate for M- 
dwarfs is consistent with the detection rate 9-23% for cold 
disks around A- to K-type main-sequence stars with ages 10 
- 180 Myr measured by Wyatt, Dent & Greaves (2003) and 
Najita and Williams (2005) in two sub-mm surveys. The con- 
sistency between these detection rates is the first indication that 
a population of cold disks is equally prevalent across stellar 
spectral types. The highest disk mass (28 ± 5 lunar masses for 
GJ842.2) among the M-dwarfs is comparable to the highest 
disk mass among the A- to K-type main sequence stars (Table 4 
and Fig. 3). We plan additional observations of a larger sample 
of M-dwarfs to ascertain this detection rate and to better de- 
termine the dust mass upper envelope of field M-dwarfs. The 
dust mass for the Kuiper Belt is ~ 3 10"^ lunar mass (Landgraf 
et al., 2002, and Moro-Martin & Malhotra, 2003), and, clearly. 



deeper observations are required to establish the statistics of 
debris disks around mature stars over a large range of ages. 
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